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Abstract. We investigated the impact of mineral dust par- 
ticles on clouds, radiation and atmospheric state during 
a strong Saharan dust event over Europe in May 2008, apply- 
ing a comprehensive online-coupled regional model frame- 
work that explicitly treats particle microphysics and chemical 
composition. Sophisticated parameterizations for aerosol ac- 
tivation and ice nucleation, together with two-moment cloud 
microphysics are used to calculate the interaction of the dif- 
ferent particles with clouds depending on their physical and 
chemical properties. 

The impact of dust on cloud droplet number concentra- 
tion was found to be low, with just a slight increase in cloud 
droplet number concentration for both uncoated and coated 
dust. For temperatures lower than the level of homogeneous 
freezing, no significant impact of dust on the number and 
mass concentration of ice crystals was found, though the con- 
centration of frozen dust particles reached up to 100 K 1 dur- 
ing the ice nucleation events. Mineral dust particles were 
found to have the largest impact on clouds in a temperature 
range between freezing level and the level of homogeneous 
freezing, where they determined the number concentration of 
ice crystals due to efficient heterogeneous freezing of the dust 
particles and modified the glaciation of mixed phase clouds. 

Our simulations show that during the dust events, ice crys- 
tals concentrations were increased twofold in this tempera- 
ture range (compared to if dust interactions are neglected). 


This had a significant impact on the cloud optical proper- 
ties, causing a reduction in the incoming short-wave radia- 
tion at the surface up to —75 W m” 7 . Including the direct 
interaction of dust with radiation caused an additional reduc- 
tion in the incoming short-wave radiation by 40 to 80 W m 2 , 
and the incoming long-wave radiation at the surface was in- 
creased significantly in the order of +10 W m -2 . 

The strong radiative forcings associated with dust caused 
a reduction in surface temperature in the order of —0.2 to 
—0.5 K for most parts of France, Germany, and Italy dur- 
ing the dust event. The maximum difference in surface tem- 
perature was found in the East of France, the Benelux, and 
Western Germany with up to — 1 K. This magnitude of tem- 
perature change was sufficient to explain a systematic bias 
in numerical weather forecasts during the period of the dust 
event. 


1 Introduction 

Aerosol particles are an important part of the atmosphere. 

They directly modify the planetary radiation budget by scat- 

tering and absorption of long and shortwave radiation, and 

they affect the properties of clouds. Depending on their 

size and chemical composition, aerosol particles can act as 

cloud condensation nuclei, CCN, and ice nuclei, IN, and 
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profoundly impacting the cloud microphysical processes and 
optical properties, hence the hydrological cycle and climate. 

The complexity of the interactions of aerosol particles with 
radiation and the hydrological cycle render the most uncer- 
tain factors in climate studies and weather prediction. 

One of the major contributors to the atmospheric aerosol 
mass is mineral dust (D’Almeida, 1986). Even though 
the emission sources of mineral dust in the atmosphere are 
mainly desert regions, dust particles can be transported over 
long distances (e.g. Swap et al., 1996; Goudie and Middle- 
ton, 2001; Tegen, 2003). Europe and the Mediterranean can 
be strongly affected by dust outbreaks originating in the Sa- 
hara Desert (Ansmann et al., 2003; Wiacek et al., 2010). Dust 
particles modify the atmospheric radiation balance (Cavazos 
et al., 2009), which can cause a cooling or heating of the 
surface depending on several factors like the surface albedo 
and the elevation of the dust layer (Cavazos et al., 2009; Yin 
and Chen, 2007; Stanelle et al., 2010). Including the inter- 
action of mineral dust with radiation in model systems used 
for numerical weather prediction has the potential to improve 
the quality of the forecast. For example in North Africa, al- 
though the background levels may be well represented by 
climatologies, there are strong dust storms affecting the ra- 
diation balance which are significantly underestimated by 
climatologies and in regions like Europe and the Mediter- 
ranean where the dust concentrations are dominated by spo- 
radic Saharan dust events (Perez et al., 2006; Chaboureau 
etal., 2011). 

Mineral dust particles are efficient IN (e.g. Field et al., 
2006; Mohler et al., 2006). Since homogeneous freezing 
of deliquescent aerosol is very unlikely for temperatures 
higher than — 35 °C, the ice nucleation in the troposphere 
strongly depends on the heterogeneous freezing involving IN 
(Cantrell and Heymsfield, 2005). In an eleven year lidar- 
derived cloud data set, Seifert et al. (2010) found a 30% 
higher amount of ice containing clouds over Europe with 
cloud top temperatures from —10 to — 20 °C in air masses 
that contained Saharan dust. Analyzing a Saharan dust out- 
break over the Atlantic in 2004, Min and Li (2010) showed 
that ice clouds tend to form at higher temperatures when dust 
particles are present. Global modelling studies demonstrated 
the important contribution of mineral dust to IN (e.g. Hoose 
et al., 2008; Lohmann and Diehl, 2006). The impact of min- 
eral dust on the formation of cirrus clouds for temperatures 
lower than — 35 °C depends on the competition of the het- 
erogeneous freezing of mineral dust with the homogeneous 
freezing of liquid aerosol particles. For conditions with low 
updraft velocities, the heterogeneous freezing of dust parti- 
cles is able to inhibit the homogeneous freezing and cause 
a decrease in ice crystal number concentration of the affected 
cirrus clouds (Barahona et al., 2010b). 

Based on satellite remote sensing, Min et al. (2009) 
showed that mineral dust significantly increases ice crystal 
number concentration and cloud droplet concentration with 
the consequence of a decreasing precipitation rate. Simi- 


lar results were found by Rosenfeld et al. (2001) with a de- 
crease in cloud droplet radius in clouds affected by dust and 
a suppression of light precipitation. In contrast, Pradelle 
and Cautenet (2002) found a decrease in cloud albedo in 
the presence of mineral dust because of an increase in cloud 
droplet radius. These conflicting results indicate that the im- 
pact of mineral dust on cloud droplet number depends on the 
individual atmospheric conditions, e.g. background aerosol 
and cloud type. A case study performed for the Eastern 
Mediterranean by Solomos et al. (2011) showed an increase 
in droplet number and in cloud top height with an increas- 
ing concentration of mineral dust particles, but also a strong 
dependency of the dust impact on precipitation on the pre- 
scribed aerosol background conditions. 

The impact of dust on the cloud droplet concentration in 
most studies was hypothesized to be caused by hygroscopic 
coatings, which make them important CCN due to their large 
size. With the exception of dust from dry lake beds, freshly 
emitted dust contains little soluble material (Kumar et al., 
2011). Processing of dust particles in clouds (Levin et al., 

1996) and the condensation of acids on the particles (Sul- 
livan et al., 2007) are prime mechanisms thought of creat- 
ing hygroscopic coatings on dust. Uncoated dust was long 
thought to be a potential CCN (e.g. Pruppacher and Klett, 

1997) , a theory sufficiently constrained by observations for 
use in models was lacking. Kumar et al. (201 1) showed that 
through the process of water adsorption, uncoated dust from 
major global source regions exhibits considerable and consis- 
tent CCN activity. In a recent global modeling study, Karydis 
et al. (2011) pointed out that even uncoated mineral dust can 
contribute significantly to cloud droplet number close to the 
source regions and that hygroscopic coating of dust particles 
can decrease the cloud droplet number in regions which are 
affected by anthropogenic emissions due to a decrease in wa- 
ter supersaturation. 

This study focuses on a major dust event that occurred 
in May 2008. Its origin was in the Sahara and from there 
mineral dust particles were transported over the Western 
Mediterranean, covering large areas of Western Europe. Dur- 
ing the episode, high aerosol concentrations were observed 
throughout Europe; IN concentrations significantly increased 
(compared to pre-event levels) at Kleiner Feldberg, Germany 
(Klein et al., 2010). During this time, traditional weather 
forecast models exhibited poor prediction skill. The Ger- 
man national meteorological service (Deutscher Wetterdi- 
enst, DWD) detected a significant bias in their numerical 
weather forecast of the 2-m temperatures of + 1 .5 K when the 
dust outbreaks reached SW-Germany (U. Damrath, DWD, 
personal communication, 2010, Fig. 1). The operational 
weather forecast at DWD is performed with the model sys- 
tem COSMO (Consortium for Small-scale Modeling) (Bal- 
dauf et al., 2011). Since the impact of aerosol particles on 
atmospheric state is accounted for in COSMO by using pre- 
scribed profiles for aerosol optical properties and prescribed 
cloud droplet and ice crystal numbers, extraordinary aerosol 
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conditions like the dust event in May 2008 are not repre- 
sented in the simulations. This can potentially cause biases 
in the weather forecast. 

The impacts of dust on atmospheric state are studied with 
the regional scale online-coupled model system COSMO- 
ART (Vogel et al., 2009) that accounts for feedbacks be- 
tween chemistry, aerosols, radiation, and clouds. A two- 
moment cloud microphysics scheme (Seifert and Beheng, 
2006) is coupled together with comprehensive parameter- 
izations for aerosol activation (Kumar et al., 2009; Bara- 
hona et al., 2010a) and ice nucleation (Barahona and Nenes, 
2009b) to simulate the impact of the various aerosol particles 
on the cloud microphysics and therefore on cloud properties 
and precipitation. This framework can quantify the impact of 
dust on atmospheric state, as it considers the competition of 
dust with other aerosol particles during cloud formation and 
the feedback processes related to cloud microphysics, radia- 
tion, and precipitation. The framework is also used to thor- 
oughly study the sensitivity of predicted atmospheric state to 
the dust amount, properties (hygroscopicity), and parameter- 
ization. 


2 The model framework 

The online-coupled regional model system COSMO-ART 
was extended with two-moment cloud microphysics and 
comprehensive parameterizations for aerosol activation and 
ice crystal nucleation. COSMO-ART is based on the opera- 
tional weather forecast model COSMO (Baldauf et al., 201 1) 
of the Deutscher Wetterdienst (DWD) which was online cou- 
pled with comprehensive modules for gas phase chemistry 
and aerosol dynamics. An extended version of MADEsoot 
(Riemer et al., 2003) is used to represent the aerosol popula- 
tion with eleven overlapping lognormal modes. Five modes 
represent sub-micron particles consisting of sulphate, ammo- 
nium, nitrate, organic compounds, water, and soot in a range 
of mixing states. These modes are coupled with the gas 
phase by condensation and nucleation, and are strongly in- 
fluenced by anthropogenic emissions of gases and particles. 
Sea salt and mineral dust particles are represented by three 
modes (Table 1) in each case. The model includes a compre- 
hensive photochemical module to capture the temporal and 
spatial distributions of inorganic and semi-volatile organic 
compounds. The emission of sea salt is calculated online 
in the model as a function of lOm-windspeed and sea sur- 
face temperature (Lundgren, 2010). The emission of dust 
particles is calculated according to Vogel et al. (2006) as a 
function of friction velocity, soil water content, and surface 
roughness by combining the parameterization of the saltation 
processes after (Alfaro and Gomes, 2001) with the treshold 
friction velocity of Lu and Shao (1999). The treatment of 
the soil parameters, which are crucial for the calculation of 
the dust emissions, are described by Stanelle et al. (2010). 
A detailed description of the treatment of gases, aerosols, and 


their emissions is given by Vogel et al. (2009) and an eval- 
uation of the simulated gas phase and aerosol characteristics 
for Europe is found in Knote et al. (201 1). 

2.1 Cloud microphysics 

COSMO includes an efficient bulk cloud microphysics 
scheme, designed for operational weather forecast (Doms 
et al., 2005). This scheme cannot treat aerosol-cloud inter- 
actions, because only one moment of the size distribution is 
calculated. Bangert et al. (2011) extended the cloud scheme 
to a two-moment representation of the cloud droplet size dis- 
tribution to consider the impact of aerosol particles on warm 
cloud microphysics. To further improve the clouds micro- 
physics and to enable the simulation of aerosol impacts on 
ice clouds a comprehensive two-moment bulk microphysi- 
cal scheme (Seifert and Beheng, 2006; Seifert et al., 2006; 
Noppel et al., 2006; Blahak, 2008) was introduced, which 
replaces the old schemes. This scheme distinguishes six 
hydrometeor categories (cloud drops, cloud ice, rain, snow, 
graupel, and hail) and represents each particle type by its re- 
spective number and mass concentration. A (generalized) 
gamma size distribution is used for each hydrometeor class, 
where the so-called shape parameters are held constant dur- 
ing the simulation. For the warm clouds, the scheme con- 
siders autoconversion of cloud droplets to rain, accretion of 
cloud droplets by rain drops, self-collection of cloud and 
rain droplets, break-up of rain drops, and evaporation of rain 
drops. Condensational growth of cloud droplets is calculated 
with a saturation adjustment technique. For the cold clouds, 
homogeneous and heterogeneous ice nucleation, diffusional 
growth of ice crystals, freezing of cloud and rain droplets, 
aggregation, self-collection, riming, conversion to graupel, 
melting, sublimation, shedding, and Hallett-Mossop ice mul- 
tiplication are considered. The freezing of cloud and rain 
drops is calculated with a classical statistical approach based 
on an empirical relation for the freezing probability as a func- 
tion of temperature (Seifert and Beheng, 2006). Note that 
the freezing of cloud and rain droplets is independent of the 
simulated IN concentrations. This may introduce an under- 
estimation of the glaciation of mixed phase clouds for high 
IN concentrations, but can be justified by a small fraction 
of cloud droplets originating from activated dust particles in 
our simulations. A detailed description of the cloud micro- 
physical processes is given in Seifert and Beheng (2006), and 
a statistical analysis of the aerosol-cloud interaction for three 
summer seasons using the microphysics scheme is presented 
in Seifert et al. (201 1). 

Although COSMO is a non-hydrostatic model and there- 
fore permits the simulation of convection, a parameteriza- 
tion of subgrid-scale deep convection has to be used for 
coarse model resolutions. We use a modified Tiedtke scheme 
(Tiedtke, 1989) for horizontal resolutions coarser than 7 km, 
which considers the convective transport of gases and aerosol 
particles. 
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Table 1. Overview of the individual aerosol modes, their chemical composition (first block), the involved cloud interaction processes (second 
block), aerosol dynamical processes (third block), and their initial diameter and standard deviation (fourth block). The model framework 
includes secondary aerosol and soot modes (if, ic, jf, jc, and s), sea salt modes (sa, sb, and sc) and dust modes (da, db. and dc). 



if 

ic 

jf 

jc 

s 
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sb 
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Activation (Kohler) 
Activation (FFH) 
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• 
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Activation (Unified) 









• 

• 

• 

Homogeneous freezing 
Heterogeneous freezing 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

Coagulation 

• 

• 

• 

• 

• 







Condensation 

• 

• 

• 

• 

• 












Initial diameter in 

jim 






0.01 

0.08 

0.07 

0.08 0.08 

0.2 

i 

12 

1.5 

6.7 

14.2 






Standard deviation 






1.7 

2 

1.7 

2 

1.4 

1.9 

2 

1.7 

1.7 

1.6 

1.5 


Table 2. Overview of the setup for the nested simulations. The 
interaction involving aerosol include anthropogenic (if. ic, jf, jc, 
and s) and sea salt (sa, sb, and sc) modes, whereas dust refers to the 
dust modes (da, db, and dc) only. 


Interactions 

REF 

R 

C 

c* 

CR 

Cloud-radiation 

• 

• 

• 

• 

• 

Aerosol-radiation 

• 

• 

• 

• 

• 

Aerosol-water clouds 

• 

• 

• 

• 

• 

Aerosol-ice clouds 

• 

• 

• 

• 

• 

Dust-water clouds (FFH theory) 



• 


• 

Dust-water clouds (Unified theory) 




• 


Dust-ice clouds 



• 

• 

• 

Dust-radiation 


• 



• 


2.2 Aerosol activation 

The activation of an aerosol particle to a cloud droplet de- 
pends on its ability to remain in a stable equilibrium with the 
ambient water vapor. The water vapor saturation ratio, S e q , 
of an aerosol particle in equilibrium with surrounding water 
vapor can be expressed as 

4crA/ w 

Seq — RT P™ D , (1) 

where cr is the surface tension at the particle-gas interface, 
M w is the molar mass of water, R is the universal gas con- 
stant, T is the temperature, p w is the density of water, and D 
is the equivalent particle diameter. The exponential in Eq. (1) 
is commonly referred to as the curvature or Kelvin effect. a w 
is the activity of the water and depends on the physiochem- 
ical properties of the particle. In case of a soluble aerosol 


particle « w can be expressed by the effective mole fraction of 
water in the solution droplet, which results in the well known 
Kohler equations. For insoluble particles (e.g. dust) this clas- 
sical theory cannot be applied, because .S' C q is now affected by 
another physical process namely the process of physisorption 
of water vapor on the particle surface. In this case cn w can be 
determined by FHH (Frenkel, Halsey, and Hill) adsorption 
theory (e.g. Fowell et al., 2004) 


— Afhh© - ®FHH 

a w — e FHH 


( 2 ) 


where Afhh. Bfhh are empirical constants, and © is the 
surface coverage (defined as the number of adsorbed lay- 
ers of water). Afhh characterizes interactions of adsorbed 
molecules with the aerosol surface and adjacent adsorbate 
molecules (i.e. those in the first monolayer). Bfhh char- 
acterizes the attraction between the aerosol surface and the 
adsorbate in subsequent layers; the smaller the value of 
■Bfhh. the greater the distance at which the attractive forces 
act (Soijamaa and Faaksonen, 2007). Afhh and Bfhh are 
compound-specific and determined experimentally. Both 
theories (Kohler and FHH) show a similar behavior for S e q, 
with a characteristic maximum (the critical supersaturation 
s c ). For ambient supersaturations greater than s c the particle 
can get activated and grow to the size of a cloud droplet. 

In this study we use a comprehensive activation parameter- 
ization (Kumar et al., 2009; Barahona et al., 2010a) based on 
a cloud parcel framework, in which a parcel of air containing 
an external mixture of Kohler and FHH particles is lifted. At 
first the s c distribution is calculated for each aerosol. Then 
the supersaturation equation of the parcel is solved using the 
bisection method to determine the maximum supersaturation, 
■Smax. and the number of activated particles, N*. As input 
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data the simulated cloud-scale dynamics (updraft, tempera- 
ture, pressure) and aerosol properties (size distributions and 
chemical compositions of all aerosol modes) are used. 

To account for hygroscopic coating of e.g. dust particles, 
Kumar et al. (2011) developed a unified activation theory by 
merging Kohler and FHH theory in their parameterization. 
For coated dust particles a w is represented by 


0! w 




A F hh© _Sfhh 


(3) 


where x w y w is the mole fraction and activity coefficient of 
water in the droplet and represents the solute effects on water 
activity. 

Based on the procedure of Bangert et al. (201 1) the activa- 
tion rate, ACT, of aerosol particles in the model is calculated 
in different ways, replacing the activation rate in Eq. (17) of 
Seifert and Beheng (2006). For a newly formed cloud the 
parameterization is directly applied and N*/At equals the 
activation rate, where A t is the time step used. In case of an 
already existing cloud the activation rate at the cloud base is 
calculated on the basis of advection and turbulent diffusion 
of particles into the cloud base 

d d d 

ACT = - — A : (Smax)w +-K—N* (W) , (4) 

OZ oz oz 

where w is the grid-scale updraft, z is the height above sea 
level, and K is the turbulent diffusion coefficient. 

In-cloud activation is calculated, if the simulated grid- 
scale supersaturation increases again above the cloud base, 
e.g. from strong updrafts. In this case, the growth of the ex- 
isting cloud droplets is considered by assuming they act as 
giant CCN that deplete supersaturation by the approach of 
Barahona et al. (2010a). The corresponding activation rate is 
calculated as A*/ At. 


2.3 Ice nucleation 


The nucleated number concentration of ice crystals, A*, is 
calculated using the parameterization of Barahona and Nenes 
(2008, 2009a,b), which is based on the framework of an as- 
cending Lagrangian air parcel. Competition between homo- 
geneous and heterogeneous freezing is explicitly considered 
in the calculation of the ice supersaturation, sp In doing 
so the dependency of A* on the conditions of cloud forma- 
tion (i.e. T , p), updraft velocity, deposition coefficient, and 
soluble and insoluble aerosol concentrations is explicitly re- 
solved. A* is given by 


jy* Ahom T Afid ( .S> Km i ) ; Ahet (.S'hom ) < A| nn 

* AhetOsj ,max ) ? -^hetCs'h om) A A| ml 

where .v M11 ax is the maximum supersaturation that develops 
in the cirrus. Shorn is the homogeneous freezing threshold 
(Koop et al., 2000), and ;V| l0m and Ahet are the number 
of ice crystals forming from homogeneous and heteroge- 
neous freezing, respectively. Ai; m is the IN concentration 



Fig. 1. Predicted as a function of forecast hour compared 
with measurements for the area of SW-Germany (see area in 
Fig. 2). Each day a forecast for three days was simulated; afterwards 
the results and the corresponding observations were averaged as a 
function of forecast hour to identify systematic differences between 
forecast and observation (U. Damrath, DWD, personal communica- 
tion, 2010). 


that completely inhibits homogeneous nucleation and sets the 
limit between combined heterogeneous-homogeneous freez- 
ing and pure heterogeneous freezing only. Ai, m is explicitly 
calculated as a function of cloud formation conditions and is 
related to the slope of the IN spectrum, Ahet(Vi) (Barahona 
and Nenes, 2009a). 

The heterogeneous freezing spectrum of Phillips et al. 
(2008) is used for AhetCsi), which is derived from several 
field campaign data sets. It provides the contribution of in- 
dividual aerosol species (dust, black carbon, and organics) 
and freezing mechanisms (e.g. immersion and deposition) to 
Ahet(si)- The contribution of the individual dust modes to 
Ahet(si) is scaled to their surface area distribution (Barahona 
and Nenes, 2009a). Niemand et al. (2012) compared differ- 
ent parameterizations of immersion freezing of dust. They 
showed that applying the heterogeneous freezing spectrum of 
Phillips et al. (2008) on dust concentrations simulated with 
COSMO-ART, the measured IN concentrations at Kleiner 
Feldberg during the dust event were reproduced. 

The homogeneous contribution to A* is given by 


Ahom — 


? 

1 

1 

1 


(Zz2A)l 

Ao 

l +e \ 7 > 


fc < 0.6 


fc > 0.6 


(6) 


where Ao is the number concentration of the supercooled liq- 
uid droplet population and f c is the droplet freezing fraction 
for cirms clouds formed in situ (Barahona and Nenes, 2008). 
Ao equals the sum of the number concentration of the solu- 
ble aerosol modes (Table 1) which are assumed to deliquesce 
during ice cloud formation. 
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The nucleation rate of ice crystals, NUC, is calculated fol- 
lowing Seifert and Beheng (2006) Eq. (34) by 

/N*-Ni \ 

NUC = max ( 1 ,0j , (7) 

where /V; is the number concentration of ice crystals before 
the nucleation pulse. 

2.4 Subgrid-scale updrafts 

Though regional models are able to resolve individual cloud 
systems, they cannot explicitly capture the updraft velocities 
which control nucleation of ice and droplets. Therefore pa- 
rameterizations of the subgrid-scale vertical velocity must be 
applied to address this issue. 

Sub-grid scale vertical velocities, w', are described with 
a Gaussian probability distribution function, P w (w'). The 
mean of P w (xv') is set equal to the grid scale updraft, ui, 
and the standard deviation, o w , is calculated as the square 
root of the turbulent kinetic energy (Morales and Nenes, 
2010). A weighted mean of the activated particles/nucleated 
ice crystals, N*, is calculated by numerically calculating the 
integral 

_ N*[ Sm .^(w')]P w W)dw' 

x / 0 °° P w {w')dw' 

2.5 Radiation 


The radiative fluxes are calculated with the GRAALS (Rit- 
ter and Geleyn, 1992) radiation scheme of COSMO for eight 
spectral bands. To consider the impact of varying aerosol 
and clouds the necessary optical properties, as they are the 
extinction coefficient, single scattering albedo, and asymme- 
try factor, have to be calculated. 

The optical properties of the aerosol particles are calcu- 
lated as a function of the size distributions, the chemical 
composition, as well as the soot and water content with the 
parameterizations of Vogel et al. (2009), Lundgren (2010), 
and Stanelle et al. (2010). 

The optical properties used for the dust particles are cal- 
culated based on measurements of the shortwave refractive 
indexes during the SAMUM (Saharan Mineral dUst experi- 
Ment) campaign and the longwave refractive indexes used in 
Helmert et al. (2007), where a mineralogical composition of 
98 % Kaolinite and 2 % Hematite is assumed. A priori Mie 
calculations for the three different dust modes were carried 
out to derive the individual mass specific extinction coeffi- 
cient, single-scattering albedo and asymmetry factor for the 
eight spectral bands of GRAALS. The actual optical proper- 
ties are calculated online for each dust mode based on the a 
priori derived values and the simulated dust mass concentra- 
tion. A detailed description of the parameterization is given 
by Stanelle et al. (2010). 

To consider the varying droplet and ice crystal size in the 
radiation scheme the parameterizations of Hu and Stamnes 
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Fig. 2. Terrain height within the model domains DO and D 1 . 


(1993) and Edwards et al. (2007) are applied to calculate the 
optical properties of the cloud droplets and ice crystals. Here 
the effective radii of clouds droplets, r c , and ice crystals, r\, 
are calculated as the ratio of the third to second moment of 
the respective size distribution. The adaption of the parame- 
terizations for the eight spectral bands of GRAALS is carried 
out following Zubler et al. (2011). Precipitating hydrome- 
teors of snow, graupel, and rain, are not considered in the 
radiation calculations, which may introduce a bias in the cal- 
culated radiation fields (Waliser et al., 201 1). 


3 Simulation setup 

A nesting approach is used for the simulations. At first 
a simulation for a domain DO covering Northwest Africa and 
Western Europe was performed. This domain enables the 
explicit simulation of the dust emissions over Africa and the 
long-range transport of the dust particles to Europe (Fig. 2). 
All gases and aerosol particles are simulated (see Sect. 2), but 
are not allowed to interact with clouds and radiation. Calcu- 
lations of the cloud microphysical processes are carried out 
with the one-moment cloud scheme, which includes cloud 
water, rain, cloud ice, snow, and graupel. The grid size is 
0.25° (~28 km) in both horizontal directions with 40 non- 
equidistant levels in the vertical up to an height of 20 km. Ini- 
tial and boundary conditions for meteorological variables are 
obtained from the IFS ECMWF analysis; for gas-phase com- 
pounds, MOZART-4 results (Emmons et al., 2010) are used. 
Clear air conditions are assumed for the initial aerosol and 
for aerosol outside of the domain. The simulation starts at 22 
May 2008, 00:00 UTC and ends 30 May 2008, 00:00 UTC. 

For the analysis of the aerosol-cloud interaction three sim- 
ulations were carried out on the nested grid for domain D1 
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(Fig. 2) covering regions affected by the dust event (France, 
Germany, and the Mediterranean). Cloud microphysical 
processes were simulated with the two-moment scheme 
(Sect. 2.1) and full interaction with the aerosol particles were 
allowed. The horizontal grid size is 0.0625°(«i7 km) in both 
directions with 40 non-equidistant levels in the vertical up 
to an height of 20 km and the simulated period is 22 May, 
00:00 UTC till 30 May, 00:00 UTC. The boundary and ini- 
tial conditions for all prognostic variables are taken from the 
simulation of the coarse grid for domain DO, which includes 
the concentration of dust transported from North Africa. The 
first four days (22 May till 25 May 2008) of the simulation 
are used to build up the aerosol concentration based on di- 
rect emissions and secondary formation. The analysis of the 
results focuses on the dust event (26 May till 29 May 2008). 
The individual nested simulations differ in the setup used for 
the dust-atmosphere interaction. In simulation REF the in- 
teraction of dust particles with clouds and radiation was not 
allowed. To investigate the impact of the dust particles due 
to their interaction with clouds we carried out simulation C 
which includes the impact of the dust particles on cloud for- 
mation via activation of dust particles and heterogeneous nu- 
cleation of cloud ice on the dust particles. In simulation C 
dust particles are assumed to not have a hygroscopic coating. 
Therefore the FHH adsorption theory is used for the activa- 
tion of the dust particles using Afhh = 2.5, and Rfhh = 1.2 
(Kumar et al., 2011). Aged dust particles can potentially 
have hygroscopic coatings which decreases the supersatura- 
tion needed for their activation. The aging of dust particles 
by condensation or chemical reactions is not treated explic- 
itly in COSMO-ART so far. For this reason, simulation C* 
is carried out in which we assume that the dust particles are 
coated by 10 % ammonium sulfate (in volume). In this case 
the activation of the dust particles is calculated with the uni- 
fied approach of Kumar et al. (201 1). In simulation CR, the 
direct interaction of the dust particles with long and short- 
wave radiation is considered additionally to their impact on 
cloud formation (with the conditions of simulation C). Ta- 
ble 2 gives an overview of the individual simulations and 
their abbreviations. 

In all simulations a detailed data set for anthropogenic 
emissions, which includes traffic emissions, emissions by 
large point sources, and area sources such as households and 
industrial areas, is used. The emission data has an hourly 
resolution and includes the emission of SO2, CO, NO x , NH3, 
32 individual classes of VOC, and the direct emission of ele- 
mental carbon particles (Pregger et al., 2007). 

4 Results 

During the 26 to 30 May 2008 Europe was affected by 
a strong Saharan dust event. A persistent trough reaching 
from Iceland to Morocco dominated the synoptic situation 
during these days. Due to steady southerly winds over 
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Fig. 3. Measured aerosol mass concentration, PMjo, at Hohenpeis- 
senberg (H. Flentje, DWD, personal communication, 2010) (black 
solid line) and simulated PM 10 of the coarse simulation for do- 
main DO (shadings) and the high resolution simulation for domain 
D1 (gray solid line, simulation CR). Red shading represents the 
total aerosol mass concentration of dust particles (mode da and 
db) and blue shading is the total aerosol mass concentration of all 
other modes (including anthropogenic aerosol particles and sea salt. 
Please note that the vertical axis is in logip-scale to take the large 
variation of PMjq during the dust event into account. 

the Mediterranean dust particles were transported efficiently 
from the emission sources in the Sahara to Central Eu- 
rope. Figure 3 shows the measured and simulated aerosol 
mass concentration of particles with a diameter below 10 pm 
(PM10) at Hohenpeissenberg, Germany (H. Flentje, DWD, 
personal communication, 2010). The measurements show 
a strong increase of the aerosol mass starting at 26 May and 
lasting till 30 May. The simulation attribute this to the arrival 
of mineral dust transported from the Sahara. The model re- 
produces the measured PM10 concentrations reasonably well 
during the dust event and captures the existing aerosol back- 
ground before the arrival of the dust which is dominated by 
particles from anthropogenic sources. Upon the arrival of the 
dust, the simulated PM10 concentration of anthropogenic and 
biogenic particles decrease from 6 pg m -3 to 1 pg m ’ and 
increase again after the dust event. This indicates that the 
aerosol population during a Saharan dust event over Europe 
differs strongly from the simple assumption of adding addi- 
tional dust particles to a typical (continental) aerosol popu- 
lation. Therefore aerosol particles have to be simulated ex- 
plicitly on the basis of anthropogenic and natural emissions 
in studies, which address the interaction of dust events with 
the atmosphere. 

A comparison of the distribution of simulated and mea- 
sured PM10 concentrations (EEA, 2011) is shown in Fig. 4 
for the main days of the dust event. The distribution of 
the dust plume (characterized by PM10 concentrations higher 
than 20pgm 3 ) is captured well in the simulations. A strong 
gradient in PM 1 0 between France and Germany is found in 
both the simulation and the measurements. On the 29 May, 
the maximum of the dust plume is shifted slightly to the 
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Fig. 4. Simulated average PM[q concentrations in the lowest model 
layer (case CR, shadings) in comparison to AIRBASE surface 
PM10 measurements (dots) (EEA, 2011) for the individual days 
during the dust event. 


North West in the simulation in comparison to the observed 
distribution. In the northern part of of the model domain and 
in the Po Valley the simulated PMio concentrations are sys- 
tematically lower, which is most likely caused by uncertain- 
ties in the anthropogenic emission inventory and unresolved 
local effects. 

At the Atlas mountains the transported dust particles were 
lifted to the upper troposphere. This ascending air flow on the 
one hand enabled the long range transport of the dust parti- 
cles to Central Europe and on the other hand caused condi- 
tions favorable of heterogeneous and homogeneous ice cloud 
formation. Additionally, warm and mixed phase clouds were 
formed over the whole domain. 

Figure 5 presents satellite retrievals of MODIS (Menzel 
et al., 2008) and calculated cloud top temperatures on ba- 
sis of simulation C during the morning hours. The cloud 
top temperature at the individual grid points are calculated as 
the temperature of the highest model layer in which the grid 
scale cloud optical thickness exceeds 0.05. Even though both 
cloud top temperatures are derived in very different ways, we 
make a coarse comparison to verify if the main features of the 
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cloud distribution are represented in the simulation. On the 
morning of the 26 May 2008 the high clouds in the eastern 
part of the model domain are well captured, whereas the high 
clouds in the South Western part are not found in the satellite 
data. On 27 May 2008, the north western part of the model 
domain is covered by high clouds which is represented in the 
simulation. The distribution of high clouds on 28 May 2008 
shows an eye-catching pattern over the Mediterranean. The 
high clouds over Eastern France and North Italy on the 29 
May 2008 are also represented in the simulation. The main 
differences on all days are in warm low level clouds which 
are detected by the MODIS retrieval e.g. over the Mediter- 
ranean but are not represented in the simulations. This dif- 
ference is partly due to a misinterpretation of the MODIS re- 
trieval algorithm which classifies the optical thick dust layers 
over water as clouds (R. Frey, CIMSS/SSEC/UW-Madison, 
personal communication, 2011). Differences in the absolute 
value of cloud top temperatures at very low temperatures 
most likely arise from the coarse vertical resolution of the 
model at these heights. 

4.1 Dust plume 

The main dust plume affecting Europe originated in Algeria 
where steady moderate winds caused an effective emission 
of mineral dust into the atmosphere. Due to south easterly 
winds, the dust particles were transported to the North where 
they were lifted into the middle to upper troposphere at the 
Atlas Mountains. 

Figure 6 presents the simulated daily mean dust number 
concentrations, Ad, for the nested model domain during the 
period of the dust event. Ad is calculated for three differ- 
ent temperature ranges, each of them favorable for a differ- 
ent process of aerosol cloud interaction. In the lowest layer, 
which reaches from the ground up to the freezing level, liq- 
uid clouds can form, with Ad being highest, up to more than 
50 cm -3 inside the plume over the Mediterranean. From 27 
May till 29 May 2008 wide areas of Eastern France and West- 
ern Germany were affected by dust concentrations above 
25 cm -3 . 

In the middle layer, which reaches from the freezing level 
up to the level of homogeneous freezing ( T ~ 235 K) and 
hence is favorable for heterogeneous freezing, dust number 
concentrations are in the order of 10-50 cm -3 with maxima 
over the eastern part of France. 

The upper layer, which reaches from the level of homo- 
geneous freezing up to 15 km height and hence is favorable 
of combined homogeneous and heterogeneous freezing, has 
mean dust concentrations in the order of 1 cm -3 over Europe. 

4.2 Dust impact on cloud droplet and ice crystal 
number concentration 

Owing to their size and surface properties dust particles 
can act as CCN. High dust concentrations can therefore 
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Fig. 5. Composites of cloud top temperatures measured with MODIS (top row) and simulated cloud top temperature (case C, bottom row). 


modify the cloud droplet size distribution and consequently 
the microphysical processes in the warm cloud phase. Dust 
particles are found to be very efficient ice crystal nuclei in 
several lab and field studies (Field et al., 2006; Mohler et al., 
2006). Ice nucleation due to immersion freezing of dust par- 
ticles starts at temperature around 263 K with a freezing frac- 
tion on the order of 0. 1 %. At temperature around 253 K, 
the freezing fraction of dust particles was found to be in the 
order of 1 % (Phillips et al., 2008). At temperatures below 
235 K ice crystal nucleation due to homogeneous freezing of 
aerosol droplets sets in. At this temperature range hetero- 
geneous freezing of the dust particles has to compete with 
homogeneous freezing for water vapor. 

In Fig. 7 joint histograms of cloud droplet and ice crystal 
number, N c and A), and mass concentrations, m c and m\, are 
shown for different temperature ranges for simulation REF 
and C. The histograms are calculated for all grid points and 
the whole period of the dust event (26-29 May 2008). With 
the help of the joint histograms, systematic differences in the 
distribution of N c , N\, m c , and m\ between the simulations 
can be identified together with difference in the overall num- 
ber of grid points containing cloud water and ice crystals in 
the individual simulations. 


For temperatures above freezing level, the domain-wide 
joint histograms of cloud droplets in simulation REF and C 
show small differences. This is because the overall num- 
ber concentration of available CCN is dominated by sea salt 
and secondary aerosol particles. The number concentra- 
tions of activated dust particles, A'j', is mostly below 3 cm -3 
and higher number concentrations up to 20 cm -3 occurred 
only for total activated number concentrations, N *, below 
200 cm -3 . The fraction of to N* is mostly below 10 %. 

Despite the fact that the coated dust particles can be acti- 
vated at lower supersaturations than the uncoated particles, 
the domain wide joint histogram of cloud droplets in simu- 
lation C* shows only slightly higher N c in comparison with 
simulation C (Fig. 7). One reason is that even though much 
more dust particles get activated in simulation C* in com- 
parison to simulation C the number is still low in compari- 
son to the total number of activated particles for most of the 
grid points (Fig. 8). But we also want to point out here, that 
there are more cloudy grid points in simulation C* where /V? 
contributes to 10% or more of N* at values of N* below 
200 cm -3 . Hence the dust particles can have a significant 
impact on individual clouds for conditions with low number 
concentrations of other aerosol particles. 


www.atmos-chem-phys.net/12/4045/2012/ 


Atmos. Chem. Phys., 12, 4045-4063, 2012 



4054 


M. Bangert et al.: Saharan dust event impacts on cloud formation over Western Europe 










26 May 2008 235 K > T 


27 May 2008 235 K > T 


28 May 2008 235 K > T 


29 May 2008 235 K > T 


26 May 2008 273K>T>235K 27 May 2008 273 K > T > 235 K 28 May 2008 273 K > T > 235 K 29 May 2008 273 K > T > 235 K 


26 May 2008 T > 273 K 


27 May 2008 T > 273 K 28 May 2008 T > 273 K 


dust number concentration in cm' 3 


o.i 


1 10 25 50 


Fig. 6. Simulated daily means of dust number concentrations for atmospheric layers below freezing level (bottom row), layers favorable of 
heterogeneous freezing (middle row), and layers favorable of homogeneous freezing (top row) for 26-29 May 2008 (columns). 


Another reason for the low impact of coated dust particles 
on A c is that the increased number of activated dust particles 
can cause lower maximum supersaturations, s max , due to the 
growth of the dust particles. This lower supersaturation can 
result in a decrease of the activated droplets from aerosols 
other than dust. This can buffer the impact of the increased 
number of activated dust particles on A c , especially in the 
case of dust particles with a hygroscopic coating (Kumar 
et al., 2009), and can be another reason of why total droplet 
number does not increase significantly in simulation C*. In 
Fig. 9 the joint histogram of the maximum supersaturation 
and the updraft during the activation process are shown for 
the cases REF, C, and C*. The impact on ,v rnax is strongest in 
case C*, with significantly lower ,v max in comparison to case 
REF Especially the number of gridpoints with .v nlax > 0.4 % 
is decreased in case C* in comparison to case REF At this 
gridpoints the concentration of aerosol particles other than 


dust is low, which explains the susceptibility of .y rnax to the 
activation of the dust particles. For uncoated dust particles 
(case C) the impact on ,v nlax is lower, because of the small 
amounts of water required for FHH activation (Kumar et al., 
2009), in comparison to the Kohler activation of hygroscopic 
particles. 

The impact of the dust particles on number and mass con- 
centration of cloud droplets and ice crystals for atmospheric 
layers with temperatures in the range of 273 to 235 K is dom- 
inated by the heterogeneous nucleation of ice crystals from 
dust particles. The joint histogram of number concentrations 
of ice crystals, A;, and mass concentrations of ice crystals, 
in j , show a systematic increase in Aj of 30 1 — 1 which corre- 
sponds to an increase of Aj in the order of 100% (Fig. 7). 
A ;Y| increase of the order 30 1~ 1 corresponds to roughly 1 % 
of /Vj which nucleate to ice crystals in the considered temper- 
ature range. Additionally, the distribution of A; is broadening 
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Fig. 7. Domain wide joint histogram of number and mass concentrations for cloud droplet (column 3 to 5) and ice crystals (column 1 to 
2) calculated for grid points with N& > 1 cm -J and different temperatures ranges for 26-29 May 2008. The top row shows results from 
simulation REF. The bottom row shows results from simulations C and C* (right figure). 
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Fig. 8. Domain wide joint histogram of activated dust particles, N 
and total number of activated aerosol particles, N£, calculated for 
grid points with N& > 1 cm -3 for 26-29 May 2008 for simulation 
C (left) and C* (right). 


from a narrow distribution of A; in case REF centered at low 
A; in the order of 1 5 1 — 1 to a wide distribution centered at 
A; in the order of 60 1 . This can be explained by the more 

variable distribution of ;V ( | in comparison to the distribution 
of soot particle IN in case REF. 

The more efficient ice nucleation in the temperature range 
between 273 to 235 K from heterogeneous freezing of the 
dust particles has an impact on the glaciation of mixed phase 
clouds. Though the freezing of cloud droplets in the realized 
simulations is independent of the simulated aerosol particles, 


the enhanced ice nucleation involving the dust particles is 
able to modify the cloud droplet distribution for tempera- 
ture below 263 K (Fig. 7). The total number of grid points 
with m c > 0 is significantly reduced in case C with respect to 
case REF. In particular, grid points with low m c and N c are 
strongly decreased in number. This is due to the Bergeron- 
Findeisen process, where ice crystals grow at the expense of 
the existing water droplets due to the difference in saturation 
ratio with respect to water and ice. The total number of grid 
points containing cloud water decreased by 61% in simu- 
lation C with respect to simulation REF in the temperature 
range between 263 to 235 K. 

For atmospheric layers above the level of homogeneous 
freezing ( T < 235 K), the joint histograms of 777; and A', are 
very similar in both simulations. At this temperature, the het- 
erogeneous freezing of the dust particles has to compete with 
the homogeneous freezing of the aerosol droplets for water 
vapor during the ice nucleation. An increasing number of 
heterogeneous freezing particles, Ahet, will at first decrease 
the maximum ice supersaturation during the nucleation event 
and therefore decrease the total number of nucleated ice crys- 
tals because homogeneous freezing is less effective. The 
value for which Ahet completely inhibits homogeneous freez- 
ing is defined as Ai; m (see Sect. 2.3). For Aj le t greater than 
(Vi; m , heterogeneous freezing is the sole contributor to ice 
crystal number concentration and will increase with increas- 
ing Ni (Barahona and Nenes, 2009a). In Fig. 10 the joint 
histogram of Ahet and Ai; m is shown for simulation REF 
and C. In the former, AV-t is mostly lower than 1 0 T 1 with 
only few grid points where Ahet exeeds Ajj m for conditions 
with Aiim lower than 60 1~ 1 . In the latter simulation, AV-t 
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Fig. 9. Joint histogram of updraft and maximum supersaturation, s max , during activation of particles calculated for grid points with 
Nd > 20 cm -3 for simulation REF (left), C (middle) and C* (right). The updraft is calculated as w + y/T KE. j max is the average of 
the corresponding updraft PDF (see Eq. 8). 


REF 

26 - 29 May 2008 235 K > T 


400 
Tj 300 
2 200 
100 


100 200 300 400 100 

Nun, in L 1 

Occurrence in 1 0 X 

1 1.5 2 2.5 3 3.5 4 

Fig. 10. Joint histogram of heterogeneously nucleated ice crystals 
and calculated for grid points with > 0. 1 cm -3 for simula- 
tion REF (left) and simulation C (right). The dashed line highlights 
the limit between combined heterogeneous-homogeneous freezing 
(below the line) and pure heterogeneous freezing only (above the 
line). 



is generally much higher due to the heterogeneous freezing 
of the dust particles, but rarely exceeds 100 1 -1 . Therefore 
(Vhet is mostly below /V| nn and is not able to significantly im- 
pact iVj. For relatively few conditions (with ;Vij rn lower than 
100 l -1 ) (Vhet is able to exceed N\ i m . 

4.3 Dust impact on cloud properties and radiation 

In the previous section we focused on the net impact of the 
dust particles on cloud droplet and ice crystal number con- 
centration for the whole model domain during the simulation 
period and found systematically higher values in ice crystal 


number concentration for atmospheric layers between freez- 
ing level and the level of homogeneous freezing in case C in 
comparison to case REF. This has numerous consequences 
on cloud properties. We already discussed the more effective 
glaciation of mixed phase clouds in simulation C, which has 
potential consequences on the vertical distribution of the la- 
tent heat release and the dynamic development of the clouds. 
The higher number concentrations in (V; (but similar values 
for m\) in both simulations causes difference in the size of the 
ice crystals and therefore has consequences on the sedimen- 
tation velocity of the ice crystals and on the optical proper- 
ties of the ice clouds. Short wave radiation is scattered more 
efficiently by smaller ice crystals, which results in optically 
thicker clouds. In Fig. 1 1 the difference in the net surface 
shortwave radiation flux, F sv/ , around noon at the 26 to 29 
May and the difference in the net surface long wave radia- 
tion flux, Fi w , around midnight at the 26 to 29 May between 
simulation C and REF is shown as a function of the verti- 
cally averaged A/j at cloud ice containing grid cells in case 
C. During noon, the distribution of the differences in F sw 
show a decrease as a function of N d reaching —75 W m -2 at 
A' ( j = 100 cm -3 . The distribution of differences in /q w show 
a slight increase as a function of A' t | . This can be explained 
by the less effective sedimentation of the smaller ice crystals 
which can cause a slight increase in the ice water path. 

The temporal evolution of the median effective radii of 
cloud droplets, r c , and ice crystals, r u calculated for domain 
Dl, in simulation C and REF is shown in Fig. 12. For r c 
slightly lower values occur in case C with maximum differ- 
ences to case REF mostly lower than 1 pm. In contrast, r\ is 
significantly lower in simulation C on the order of 10-25 pm. 
The difference in r\ scales with the time evolution of the mag- 
nitude of r\. Small r; indicate that ice clouds occur mainly at 
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Fig. 11. Joint PDFs of differences (C minus REF) in surface net 
shortwave radiation flux around noon (left) and in surface net long- 
wave radiation flux around midnight (right) together with average 
dust concentration at grid points containing cloud ice (C). The 
dashed line is the result of a linear regression fit of the data. 

greater heights, where the difference in r\ is much smaller 
due to the controlling influence of homogeneous freezing. 

Figure 13 shows the temporal evolution of domain aver- 
aged cloud properties and radiation fluxes in simulation C 
and REF. The temporal evolution of averaged vertical in- 
tegrated cloud water content, LWP, and ice water content, 
IWP, shows values of similar magnitude in simulation C and 
REF with positive and negative difference of a few percent. 
Whereas LWP shows most of the time slightly lower values 
and IWP slightly higher values in case C. Differences in LWP 
and IWP are mainly caused by differences in the dynamics 
of the cloud systems. The average in-cloud vertical velocity, 
Wcioud, shows slightly higher values in case C during the first 
three days of the dust event. This can be explained by more 
effective glaciation of mixed phase clouds in case C and the 
additional release of latent heat. On the last day w c i ou d is 
lower in simulation C which is also reflected in a lower LWP. 
In general, differences in dynamics can be driven by different 
processes, e.g. changes in radiation and consequently tem- 
perature caused by modified cloud properties, and are hard 
to relate directly with a specific individual process. 

In both simulations, the average hourly precipitation rate, 
R\ h, is almost identical. Difference in R\ h scale mostly with 
the differences in LWP and IWP, e.g. during the last two days 
of the dust event were the largest differences occur. Despite 
the almost identical R i h , the local differences can be large 
due to spatial shifts in the distribution of precipitation. 

The domain average F sv/ shows systematic differences be- 
tween simulation C and REF with lower values in simulation 
C. The maximum F sw is reduced significantly in case C by 
10 to 20 W m -2 . The largest difference in F sw occurs on 27 
May with 20 W nr 2 lower values in case C. During this day 
the dust amount is highest in the model domain. Because 
LWP and IWP are almost identical in both simulations, the 
difference in F’ sw can be directly related to the impact of the 




Fig. 12. Time series of the average effective radii of cloud droplets 
r c (top) and ice crystal r; (bottom) of domain D 1 . Solid lines are re- 
sults from simulation C and dashed lines are results from simulation 
REF. r c is calculated for gridpoints with a cloud water mass mix- 
ing ratio greater than 0.05 g kg - 1 and r; for ice crystal mass mixing 
ratios of more than 0.00 1 g kg - 1 . 

dust particles on the radiative properties of the ice clouds in 
case C (Fig. 11). The temporal evolution of the average /j w 
shows just small differences between simulation C and REF 
with a maximum of +1.5 W m“ 2 . 

The difference in the net radiation fluxes at the surface 
have an impact on the simulated temperatures. The domain 
averaged daily temperature maximum of the temperature in 
a height of 2 m above ground, T 2 m is up to 0.3 K lower in 
case C than in case REF. Tj m is lower in case C until the 
daily temperature minimum is reached in the early morning, 
where Tim in both simulations converge again. Despite the 
small difference in the domain average 72 m , the differences 
are locally much larger. This will be discussed in Sect. 4.5. 

4.4 Direct dust impact on radiation 

Until now we investigated only the impact of dust on the at- 
mospheric state due to the interaction of dust particles with 
cloud microphysical and optical properties. During the dust 
event the dust particles interact also directly with the radia- 
tion fields of the atmosphere due to scattering and absorption 
of short and longwave radiation. Results of simulation R, 
which considers only the direct interaction of dust particles 
with radiation, and simulation CR, which includes addition- 
ally to the interaction with the clouds also the direct interac- 
tion of dust particles with radiation, are shown together with 
the results of simulation C and REF in Fig. 13. The domain 
averaged dust optical depth at a wavelength of 550 nm, rjust. 
increases during the dust event up to 0.36 at the 27 May. Dur- 
ing the last two days of the dust event td ust decreases again 
to 0.1. The additional direct interaction of the dust particles 
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Fig. 13. Time series of dust optical thickness (rdust). cloud properties (u) c i ouc j, LWP, and IWP), precipitation rate (R\ h), net surface radiation 
fluxes (F sw and Fj w ), and temperature (72m) averaged for domain Dl. Solid lines are results from simulation C, dashed lines are results 
from simulation REF, solid gray lines are results from simulation CR, and dotted gray lines from simulation R. 


with radiation in case CR causes lower values of F sw of up 
to -80Wm^’ at 27 May in comparison to case REF. On 
the other days of the dust event the maximum F sw is 40 to 
60 W m -2 lower in case CR with respect to case REF. Due 
to the interaction of the dust particles with the long wave ra- 
diation Ei w is in the order of 3 to 10 W m -2 higher in case 
CR in comparison to the cases C and REF. The differences in 
F sw and /’i w between case CR and case C and the differences 
between case R and case REF are consistent. This shows that 
changes in cloud cover and cloud properties induced by the 
direct interaction of the dust with radiation are negligible in 
the presented case. The direct impact of the dust particles in 
simulation CR causes lower daily maxima of 72m when com- 
pared to simulation C. The domain averaged 72m is 0.5 K 
lower in CR in comparison to REF at the 27 May. During 
the night T 2 m converge in the different simulations and get 
even slightly higher values of 72 m in case CR in comparison 
to C due to the higher values of F| w . The temperature de- 
creases during the day in R in comparison with REF (by the 
direct radiative impact of dust) and in C (by the interaction 
of dust with the cloud properties) in comparison to REF are 
comparable in magnitude. The sum of both effects is compa- 


rable with the temperature decrease in CR in comparison to 
REF (by the combined direct radiative impact of dust and the 
interaction of dust with the cloud properties). 

Whether the direct impact on the radiation or the interac- 
tion with the clouds has the larger impact on the domain av- 
eraged 72m varies from day to day depending on the spatial 
distribution of clouds and dust particles. 

4.5 Impact on 72 m 

In Fig. 14 the impact of the dust on average 72m is depicted 
for interaction with cloud properties only (C minus REF) and 
for direct interaction with radiation together with the interac- 
tion with cloud properties (CR minus REF) for the whole 
model domain. The 72 m for case C at noon is systemat- 
ically lower over the continent by —0.2 to —IK. Isolated 
areas where 72m in case C is higher than in case REF can 
be attributed to small spatial and temporal differences in the 
distribution of the clouds. The difference in 73 m at midnight 
is lower in comparison to noon time. The 72 m at midnight 
in case C is lower by —0.1 to — 0.4 K over the continent. 
The maximum difference is in the eastern part of France. 
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Fig. 14. Mean difference in T^m at 12:00 UTC and 24:00 UTC from 
26-29 May 2008 and total mean difference in 72 m from 26-29 May 
2008 between simulation C and REF (left column) and simulation 
CR and REF (right column). 


On average over the whole dust event (26-29 May) Ti m in 
case C is systematically lower than in case REF by —0.2 to 
—0.5 K over the continent. The maximum differences in 7 2 rn 
are in the areas of the maximum average dust concentrations 
(Fig. 6). 

The average difference in 72 m from combined interaction 
of dust with clouds and radiation (difference between simu- 
lation CR and REF) shows negative values during noon time 


of up to —1.8 K in the areas of the maximum dust concentra- 
tions. During midnight, the average difference in 73 m is only 
negative in the areas with a high negative difference during 
daytime. In all other areas the difference is slightly positive 
with values up to + 1 K locally. The average difference in 
72m between the simulations CR and REF calculated for the 
whole period of the dust event shows systematically negative 
values in the order of —0.2 K for the whole continent with 
maxima up to — 1 K locally. 

As mentioned before, an analysis of the weather forecast 
results for SW-Germany showed too high '1 2 m (by up to 
+2 K) in comparison with station observations in the after- 
noon during the dust event. In the early morning hours, the 
predicted 72 m converged with the observations (Fig. 1). 

In Fig. 15 the temporal evolution of the simulated 7) m av- 
eraged over an area including SW-Germany, G, is shown for 
the cases REF, C, R, and CR. The differences in m be- 
tween C, R, and CR and REF reach their maximum values in 
the afternoon. In the early morning hours T 2 m in the different 
simulations converge again. Beside the 28 May, Tjm in case 
R is higher (+0.1 K to +0.3 K) during the night in compari- 
son to REF due to the higher incoming longwave radiation at 
the surface (Fig. 13). The differences between C and REF, 
and R and REF are of similar magnitude during day time, 
whereas the impact of dust on 72m was dominated by the in- 
teraction of dust with clouds at the 26 May and 29 May and 
by the direct interaction of dust with radiation at the 27 May 
and 28 May. The difference between CR and REF reaches 
—0.65 K to — 1 K at the afternoon, which is consistent with 
the sum of the differences between R and REF, and C and 
REF most of the time. This indicates that feedbacks between 
the direct radiation impact of dust and clouds are of minor 
importance in this case. 

The simulation results presented in this study indicate that 
the bias in the operational weather forecast can be attributed 
to the missing interaction of the dust particles with clouds 
and radiation. Although the observed bias of 72m in the 
weather forecast is two times larger than the simulated max- 
imum difference in 72m in the area Q; the daily cycles are 
comparable. The differences in the magnitude can be ex- 
plained by uncertainties in the cloud cover and in the ex- 
act position of the dust plume; e.g. the differences in 72m in 
Eastern France reach magnitudes comparable to the observed 
bias. 

Both processes (the interaction with clouds and the direct 
interaction with the radiation) contribute in a similar magni- 
tude to the simulated differences in 72m- Whether the direct 
impact on the radiation or the interaction with the clouds has 
the bigger impact on 73m varies depending on e.g. the spatial 
distribution of clouds and dust particles. 
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Fig. 15. Temporal evolution of 72 m as simulated by the different simulations (top left) and difference in 73 m of the individual simulations 
with respect to case REF for area (see Fig. 14). Additionally to the difference between CR and REF (bottom right) the sum of the 
differences between R and REF, and between C and REF is shown (dotted line, C+R). 


5 Conclusions 

We investigated the impact of mineral dust particles on 
clouds, radiation, and atmospheric state during a strong Sa- 
haran dust event over Europe in May 2008, applying a com- 
prehensive online-coupled regional model framework that 
explicitly treats particle microphysics and chemical compo- 
sition. Sophisticated parameterizations for aerosol activation 
and ice nucleation, together with two-moment cloud micro- 
physics are used to calculate the interaction of the different 
particles with clouds depending on their physical and chem- 
ical properties. 

It is shown that the model framework is able to reproduce 
the measured aerosol mass concentrations during the dust 
event reasonably well and capture the aerosol background 
before the arrival of the dust, as well as the spatial distri- 
butions of the clouds in comparison with satellite measure- 
ments. Dust particles act as CCN and IN in the atmosphere 
and can interact with clouds in several ways depending on 
the atmospheric conditions and cloud type. 

For temperatures lower than the level of homogeneous 
freezing the dust particles have to compete with the homoge- 
neous freezing of liquid aerosol particles for water vapor dur- 
ing the ice nucleation. Though the concentration of frozen 
dust particles is up to 100 1“ 1 during the ice nucleation events 
we found no significant impact on the number and mass con- 
centration of ice crystals in this temperature range. 

The impact of the dust on cloud droplet number concentra- 
tion was found to be low, with just a slight increase in cloud 
droplet number concentration for both uncoated and coated 
dust. The number of activated dust particles was found to be 
in the order of 1 to 20 cm -3 , with higher numbers in the case 
of the coated dust particles. The additional activation of dust 
particles caused lower maximum supersaturations during the 


activation process, especially in the case of the coated dust 
particles. 

Mineral dust particles are found to have the largest impact 
on clouds in a temperature range between freezing level and 
the level of homogeneous freezing, where they determine the 
number concentration of ice crystals due to efficient hetero- 
geneous freezing of the dust particles and modify the glacia- 
tion of mixed phase clouds. Our simulations show that dur- 
ing the dust events ice crystals concentrations were increased 
twofold in this temperature range in comparison to a refer- 
ence simulations which neglects the interaction of the dust 
with the atmosphere. As a consequence the liquid water frac- 
tion was reduced significantly with 60 % less grid points con- 
taining cloud water in the temperature range between 263 K 
to 235 K. 

The strong increase in ice crystal number concentration 
has an influence on the optical properties of the ice clouds 
during the time of the dust event due to a significant de- 
crease of effective radii of the ice crystals in the order of up 
to —25 pm on average. It was shown that depending on the 
number concentration of dust available for ice nucleation the 
incoming short-wave radiation at the surface was reduced up 
to —75 W rrT 2 at Ad = 100 cm -3 . Consequently, a reduction 
in surface temperature in the order of up to 1 K over Europe 
during the day time arises from aerosol-cloud interactions. In 
the morning hours, the surface temperatures converge again 
between the different simulations. On average over the whole 
period of the dust event a reduction in surface temperature in 
the order of — 0.2 K to —OAK was found in the eastern part 
of France and Western Germany. 

The simulated aerosol optical thickness of the dust parti- 
cles was found to be in the order of 0.2 to 0.5 on average 
during the dust event over Europe. This caused an additional 
reduction in the incoming short-wave radiation which was 
found to be in the order of — 80 W m -2 on 27 May. On the 
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other days of the event the reduction was found to be in the 
order of —40 to — 60 W m~ 2 . In contrast to the aerosol-cloud 
interaction only, the incoming long-wave radiation at the sur- 
face was increased significantly in the order of + 10 W m -2 . 

The direct impact on radiation together with the impact on 
the clouds caused lower surface temperatures of up to — 1 .8 K 
in the areas of the maximum dust concentrations during the 
daytime. Due to the enhanced incoming long-wave radiation 
a temperature increase was found for most parts of Europe 
during the night time. Only in the areas with the maximum 
temperature difference during the day the night temperatures 
were still below the ones simulated neglecting the interaction 
of dust with the atmosphere. 

On average the overall impact of the dust caused a reduc- 
tion in surface temperature in the order of —0.2 to —0.5 K 
for most parts of France, Germany, and Italy during the dust 
event. The maximum difference in surface temperature was 
found in the East of France, the Benelux, and Western Ger- 
many with up to — 1 K. 

The simulations suggest that the bias which was found 
in numerical weather forecast during the period of the dust 
event can be explained by the missing aerosol cloud interac- 
tion in current weather forecast simulations. Weather fore- 
cast models use climatological aerosol conditions to repre- 
sent AOT and calculate ice crystal and clouds droplet number 
concentrations. Hence they are not able to represent condi- 
tions which differ significantly from the climatological val- 
ues, which is the case during strong dust events over Europe. 
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